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SYNOPSIS 


In this work critical speeds of rotating shafts 
are obtained. Isotropic and anisotropic (dissimilar 
stiffness in two principal directions) shafts are onsidered. 
Various factors affecting the critical speeds are taken into 
account and their effects are studied. 

The relevant governing differential equations are 
derived and FEM has been used for their solution, *.G-alerk'iif's 
technique is used to minimise the residues obtained by finite 
element approximation. The element matrices are integrated 
numerically using Gauss quadrature, assembled and solved for 
frequencies and critical speeds for various boundary conditions. 

The factors affecting the critical speeds considered 
for the present study are shear deformation, rotatory inertia 
and gyroscopic moments. Shaft with dissimilar support 
conditions in two principal planes has also been studied. 

;*n attempt to estimate the errors due to the assumptions of 
massless shaft and rigid disk has been dene. The critical 
speeds for multiple span shafts are also obtained. 

The effect of internal damping on the stability 
has been studied. Both viscous form and hysteretic form of 
internal damping are considered* The effect of bearing , chare - 
cteristics (stiffness and damping) on the instability threshold 


has been studied 



CHiYPTEtl- I 


IfJTRODUCTION 


1 .1 TrlE PROBLEM IN BRIEF : 

It was recosnised early in 18th century that, 
a rotating shaft initiated vibrations of large amplitudes 
at certain speeds and often caused catastrophic failure of 
supporting structure and other parts of the system. The 
operating speeds at v/hich this phenomenon occured, came to 
be known as critical speed, 

A’ shaft always has certain amount of mass unbalance, 
which during rotation, deflects the shaft axis from its 
equillibrium position. This leads to transverse vibration 
of shaft at its natural frequencies. In general the natural 
frequencies of a rotating shaft are dependent on speed of 
rotation. The vibrations of shaft along any two mutually 
perpendicular directions in the plane of cross section, 
occur at same frequency, but with a phase difference of + 

Due to this^the shaft axis begins to process about the line 
joining the supports, with a speed equal to the anguler 
frequency of natural vibrations . The precession of shaft 
can be in the same direction as that of spin of shaft or in 
ithe opposite direction. V/hen precession is in same direction 
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as llir-'t of spin, it is iErmed as forward precession and 

when in opposite direction, it is called as baclward precession. 

The inertia force due to unbalance has the frequency 
of spin speed of shaft. Hence, v/hen the speed of precession 
of shaft becomes equal to the spin speed, resonant state is 
reached, i,e, amplitude grov/s linearly v/ith time. When this 
state arises, the corresponding spin speed of shaft is 
called as critical speed. Depending upon the direction of 
precession, it is termed as forward critical speed or 
backv/ard critical speed. 

1 ,2 LITSRATURS REVIEW ; 

The first published work related to this problem 
v/as in 1869» by Rankine [l ], The operating speeds of rotating 
machinery were below the first critical speed in those days . 
Hence, most of the early v/ork in this field is limited to the 
determination of first critical speed. However, with the 
advent of new technology in the field of rotating machinery, 
the present day operating speeds are fsr in excess of first 
criticcl speed. Consequently, the more recent literature 
treats a greater range of problems related to rotating shafts,, 
such as effect of various system parameters on critical 
speeds, determination of stable operating speed ranges, 
causes of instabilities and methods to improve the stabilities , 
Thus the study of vibratory motion becomes necessary, for ,* 
stable and safe operation of rotating machinery. 
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A very good survey of v/ork done in this field till 
1967, can be found in the monograph ^Dynamics of rotating 
shaft" by Loev;y and Piarulli ]. This monograph tries to 

illustrate the basics of rotor dynamics. It neglects the 
mass of the shaft, i .e it uses lumped mass approach. It provides 
general but qualitative information on the effects of 
V ,rious parameters on critical speeds ‘and stability. It. 
also treats coupled, lateral and torsional motion and 
balancing . 

A more involved and analytical study of the problems 
of rotor dynamics is found in the work of Tondl [3 ] . The 

book gives a good insight of several phenomena of rotor 
dynamics. The topics emphasised in his v:ork are, effect of 
internal damping, effect of unequal shaft stiffness and effect 
of oil film in Journal bearings. It deals comprehensively with 
problems concerned with the stability of rotor motion, v/ith 
the creation of self excited vibration, and_with non-linear 
resonances. But almost the whole book is based on the 
assumption of massless shaft and rigid disk. All the same, the 
book gives an excellent reviev; of the problems, supplemented 
with results of experimental investigations. Combined effect 
of torsional and flexural rigidities, combined effect of 
several factors, causing instability have been dealt in 
detail . 
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^inother very prominent contribution in this field 
is of himentberg 4 j , This worl: deals v.'ith fle>cural 
transvdi'se vibration. Major part of the book uses continuous 
approach, closed form solutions, for hinged -hinged support 
conditions, for several cases are obtained. This book shov/ed ^ 
the limitation of critical speeds in i^rw'ird precession of a 
shaft v/ith uniformly distributed mass,v,'cs a result of 
neglecting the effect of shear deformation. Effect of 
external and internal damping, effect of dissimilority of 
stiffness, effect of dissimilar fixing conditions in two planes, 
multi -bearing systems, are feu of the cases studied in his 
work. Further, both uniform stiffness and dissimilar stiffness 
shafts are analysed with and v/ithout consideration to shear 
deformation effect. 

Another recent text in this field is of Rao [5 ] 

It uses lumped mass approach. 

An excellent collection of papers on several 
problems of rotor dynamics, presented at an international 
symposium held at Denmark in 1974, can be found in the book 
edited by Niordson [6 ]. 

The application of FEM to the problems related to 
rotor dynamics was done by Ruhl and Booker [7 ] » ii^ 1972. 

They analysed stability and unbalance response of a turborotor 
system. But the, model used, included only translational 
kinetic energy and elastic bending energy i,e,. used Euler 
beam theory. Nelson [3' ] improved the model by incorporating 
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the effect of rotetory inortio and gyroacoiic action. 

Formulation is done using energy method. The effects of bearing 
stiffness and damping are also incorporated. Zorzi and 
I'Jelson [ 9 ] , incorporated the effect of internal damping 
in the finite element model of [ 8 ] • --ccording to [4 ] 

the possibility of instabilities in the post -critical speed 

) 

rtange due to internal damping, was first shown by Kimball 

and Niwkirk [ 10 ] * mechanism of internal damping is 

better represented by considering the hysteretic form, rather than 

viscous form [ 4 ] • combined model incorporating both the 

forms of damping is formulated by [ 9 ] • ^ study of 

effect of bearing stiffness and damping charecteristics on 

the instability threshold is done by them. Later in 19S5, 

Greenhill, Bickford and Nelson [ 11 ] , extended the work 

of [ 9 ^ , for a shaft with variable cross section, by using 

a conical beam finite element. 

There are a good number of authors who use transfer 
matrix method in this field. Lund [12 ]» studied the 

effect of fluid-film journal bearings and hysteretic internal 
damping, on the critical speeds and instebility threshold. 

He also incorporated the effect of bearing stiffness and 
bearing damping. But he treats the shaft as a uniform beam, 
neglecting the effect of shear deformation , rotatory inertia 
and gyroscopic moments. Bansal and Kirk [13 ]» improved this 
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■'•■(Ork by incorporcting the disk gyroscopic moments and 
shear deformation effect, Ramnlo’ishnan and Prabhu [l4 ] 
have studied the unbalance response, instability and critical 
speeds of multispan rotor systems , 

■lO-l these authors, including Rao [ 5 ^use lumped mass app- 
roach in the transfer matrix method . Lumped mass approach is 
* al\.’a3?s less accurate than the distributed mass approach. 

Moreover this method requires some skill and experience as 
there are chances of missing some modes , Satisfaction of 
boundary conditions of intermediate supports is not very 
straight forward, 

Vi’ith continuous approach in transfer matrix method 
for beam problems^ one faces some numerical problems. 

Moreover it is quite difficult (if not impossible) to obtain 

transfer matrices for rotors v/ith rotatory, shear and 

* 

gyroscopic effects , This difficulty is easily overcome in 
transfer FEM (TFEM) at the cost of little accuracy, Gupta [ 15 ] 
and Goel [16 ] .It eliminates the need of the solution 
over an element as it uses the finite element shape functions , 
Further, TFEM also has the advantage of lesser computer 
memory requirement as, the order of matrices to be handled 
remain same as in transfer matrix method , Thus the TFEM blends 
the advantages of FM and transfer matrix method. But 
obtaining the complex frequencies for damped cases does not 
seem possible. 
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uecently, Goal [ 16 ]obtained the criticcl speeds of 
uniform, isotrordc undamped shafts usin^; Transfer FEM. The 
effects of shear deformation, rotatory inertia and gyroscopic 
moments were included in his analysis. It also covers 
multi-span shafts . In his v/ork, two second order governing 
differential equ-'tions ^solved using linear polynomial, 
satisfying the continuity and completeness requirement led to 
ill- conditioning . To avoid this, these two equations v/ere 
combined into -a- fourth order differential equation and 
cubic polynomial used. Although it was not possible to 
satisfy all the geometric and natural boundary conditions 
exactly in this approach, the results \vere quite 
satisfactory . 

1 .3 AIM AW) SCOPS OF PRESENT WORK ; 

The present v/ork aims at obtaining the critical 
speeds of flexible rotating shafts and studying the effect 
of various parameters on critical speeds , The relevant 
governing equations were derived and Galarkin FEI''I has been 
used to solve these more farailier differential equations . 

The effects of gyroscopic moments, rotatory inertia 
and shear deformation are incorporated. A study of critical 
speeds, neglecting the shear deformation is also done. Two 
second order governing equation approach is used in the present 
work, except for the cases studied neglecting the effect of 

' t 

shear deformation. 
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For the.' vibr"tio:'3 analysis oi sha.fts, carrying disks 
having m.ss larger than shaft mass, the mass of shaft is 
generally neglected. An attempt has been made to estimate 
the error committed due to this assumption. For this, critical 


speeds for a range of ratios of shaft mass to disk mass 
are obtained. Further, the effect of flexibility of disk is 
also studied. 

The effect of dissimilar shaft stiffness, caused by 
the unequal area moments of inertia of shaft section in the 
two principal planes is studied. The effect of dissimilarity 
in support conditions in the two planes is also included. 


A study of instabilities caused by internal damping 
and the effect of bearing stiffness and damping properties 


on the stability ahe studied. Hysteretic, as well as viscous 
forms of damping are considered. The combined model 
incorporating both the forms, is formulated bssed on the work 
of I'lelson [9], The model is studied independently for both 
the forms of damping, Variation in instability threshold, 
with bearing stiffness and bearing damping is also studied. 



CHAPTER- II 


THEORY A^D FORMULATION 


This chapter gives the derivation of governing 
differential equations and corresponding FEM formulation 
of the various cases of shafts studied in this work. 


2,1 ISOTROPIC SHAFT WITH NO DAMPING : 


Governing differential equations for transverse 
vibrations of rotating shaft, incorporating the effects of 
shear deformation, gyroscopic moments and rotate r^^ inertia 
are derived from Hamilton's principle. 

Consider a differential element of length ds, as shown 
in Fig. (2.1). The kinetic energy of the shaft becomes [ S ]. 


Tg=/^ (x^+y^)ds+ / ~ ®y^ ® 3 c®y 

Here, m is the mass per unit length ; x and y are tbe 


translations along OX and OY axes ; 1^ is the diametral mass 
moment of inertia per unit length . 0^ and ©y are the bending 
slopes in YS and XS planes; 1^ is the polar mass moment of 
inertia per unit length and is the spin speed, assumed 
constant. Dot represents the derivative with respect to time t. 
In the above Eqn; (2,1), ‘the first term is the energy due to 
translation, second term is the energy due to rotatory inertia, 
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X 


Fig. 2.1 Free body diagram of differential 
element in two planes 
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tlie third terra is the energy due to spinning of the shaft and 
the last term gives the energy due to gyroscopic moments, 
derivation of kinetic energy is given in Appendix I , 


The strain energy of the shaft is 



( 2 . 2 ) 

where El is the flexural rigidity ; K* is transverse shear 

form factor, a constant depending upon the shape of cross 

section } G is the shear modulus ; A is the area of cross section; 

/? and /? are the angles of distortion due to shear deformation 
X y 

about X and Y axes respectively. Strokes represent the 
derivatives with respect to S. First tern in Eqn. (2.2) 
represents strain energy due to bending and the second term , 
strain energy due to shear deform*: tion. 


It is seen from the Fig, (2.1 Ithat total slopes are 


X = 


y' = 


V^y 


♦ • • • 


(2.3) 


Using Sqns. (2.3), expression for strain energy U given by 
Eqn, (2.2) becomes 

1 1 ^ ? 

U= ; 1 EI(q‘ ^+ 0' ^)ds + / ^ GA[ (x' -0 )^-t-(y» - 0 ^) ] ds 

Q ^ ^ y o 

.... (2.4) 

Thus integral to be minimised is 

I = ! (I-U) « 
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1 


1 r '2. , 2 ^ 1 

_ I tr \ 1 ' 


^ \ i m(x''+y"^)+ ^ I^^(0^+e2)+ ^ I, 

'1 


t 0 *'2 I ^ ^p C o 


- S °o ®x®y- 1 El 6.2, 

-jK'GA f.x' -9y)2+(y' -6^)^ ]1dsdt 


By Hamilton’s principle 

tp 

6I = 6 / (Tg-U)dt = 0 


.... (2.5) 


.... ( 2 . 6 ) 


Applying operator 6 to Eqn, (2.5) is quite straight 
forward, but fourth term needs little special attention, and is 
^2 ^ 


^ Q„ 9„ ©„ ds dt 

t^ 0 p ° ^ y 

1 tr 


"^2 ^ 




0^6©^ ds dt 


/ I C, 60 © 
n P O ^ y 


1 t 


ds- / 1 66^6 dt ds 

O 


'2 1 


* I ®x *®y <1® ** 

o 


t^l 


= i ®o <-®y °®x+ ®x ^®y) 

t.^ o 

.... ( 2 .. 7 ) 

Thus. applying the 6 operator to all terms in Eqn. (2.5) 
and rearranging, one gets 
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1 , 
j" S x+(K'GA(x' -0 ) ) ]5x 

o ^ 

[-^D V^p +K'G«{y'- e^) ]6e^ 

+ [-my + (k'gaCv' -9^^))' ]i5y 

+ [ -Ij, 9y-Ip 9^+(EI 9^' +k' GA{x'- 9y) ] 69y] ds dt 

- El e;, C -El 9 ; 69 

0 o 

1 1 

-K'GACx' -0 )6x -k'GA (y'-© ) 6y =0 
^ o o . 

•••• (2»8) 

The displacements 6x , 6© 6y and 6© are arbitrary, therefore 

X y 

(K'Ga (x'-©y)) - mx = 0 .... (2.9) 

(El ^y)' + K'GA (x- -9y) - Ij) 9y- Ip 0p9^ = 0 

(K'GA (y*-0jj)) - my =0 .... (2,10) 

(El ©y +K'GA(y‘ - ©^)-Ij) VIp ^0 ®y = ^ 

El ©^ 6©^ P = 0 ; k'gA (x*- © ) 6x I = 0 

.... (2.11) 

|1 I ^ 

El ©' <5©^ = 0 ; K GA (y* - © )' 6y =0 

( 2 . 12 ) 



Equations (2.9) are the governing differential equations in 
:LS plane, Eqn, (2.10) in YS plane and Eqns . (2.11) and (2.12) 

are the corresponding boundary conditions. These equations can 
also' be derived by Newton’s method [4 ] , (see Appendix II), 

Multiplying Eqns, (2.10) by i = , combining it v/ith 

Eqns. (2,9) and introducing the expressions 

2 = x+iy and ^ = 0^ + (2.13) 

complex differential equations for isotropic shaft are obtained as 

(k'ga(z' -{?())' - mz = 0 .... (2.14) 

(El )' + k' GA (z' - ) -Ij3 5^ + ilp {5 = 0 

and for isotropic boundaries the complex boundar^^ conditions 
become 

El Cf' = 0 or 0 s= 0 at X = 0 and 1 
, .... (2.15) 

K GA(z‘ - 0 ) = 0 or,z=0 at x=0 and 1 


The two equations given above can be combined into a single 
fourth order differential equation as 


El z 


^"^-It-XI + -~) z*'+ i Qq ^‘+niz + 

JJ Q p o 


ra It, ml„ Q- 

« i P ° 

y! GA 


K? GA 


2 = 0 


.... (2.16) 

which is same as that given by [ 4 ]. 

If shear deformations are negligible, the above 


equation becomes^ 
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El z" + 1 Ip 0^z“ + jfi = 0 .... (2.17) 

The second and third terms represent the rotatory inertia and 
gyroscopic actions respectively. The above equstion is same 
as that given in[ 3,A ]. 


2,1,1 Mon dlmens ionalisation : 

Ibie governing differential equations (2,1 A) and 
(2,17) are non-dimensionalised using the follov/ing parameters. 


Z = 


A = 


I * 


c S 

^ - T i 


T - ^ 

T 


Q = 


* 


Q^T 

o 


n 


. T I 


* PJ 


P^o^ 


( 2 , 18 ) 




Vvhere and are the area and area moment of inertia in 
o 0 

bending at S = 0. 

Using these parameters, Eqns, (2,14) and (2,15) 


become 

(A^ (Z' - ^) )’ - Me Re^ Z = 0 - (2.19) 

(I 0 ^ )' + 1_ A (z‘ - CJ)- In Re^ 0 +i 21 Re^ Q {5 = 0 

^ Me Re^ ^ 


Me Re^ I^ 0'= 0 

or 0=0 

(2.20) 

- 0 ) = 0 

or Z = 0 


where Me = — and Re^ 
K'G 

^0 

(2.21) 


It may be noted that now prime (;) indicates 
differentiation with respect to S and dot (,) with respect to T, 



For uniform shaft and are unity , and 
Eqns. (2.19) and (2.20) become 

(Z' - 0 y - Me Re^ Z = 0 




// 1 


5“ (Z' - (Zl)- Re^ ? + i 2Re^ Q 0 


.... ( 2 . 22 ) 


McRe 


Me Re^ 0 or 0 = 0 

(Z' - (25) = 0 or Z = 0 

Similarly Eqn, (2.17) becomes 


- Re^ & - 2i 0 z") + Z = 0 


.... (2.23) 


.... (2.24) 


2.1.2. FEM Equations : 

Equitions (2.19) are partial differential equations. 
Their solution is obtained by the finite element method , This 
method leads to simultaneous ordinary differential equations 
which can be solved in a routine way. 

Descretizing the domain (0,1) by number of finite 
elements of length h (Refer Fig, 2.2) the variation of (^ and Z 
over a typical element can be approximated as 

=[W^J{Z) ....(2.25) 


= a'+b' 0+ = Ln®* J 

where [K^ J and [N® J are the shape functions j tZ} 
are the nodal parameters. 


(ne) 


and 




Fig. 2-2 Finite element discretization 
and typical element 


1 
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Substituting Eqns. (2.25) in Eqns . (2.19), residues 
obtained as 


are 




R^ = 


(A„ ))' . MO Re^ Z<-^'> 

)'+ — - — 2 .... ( 2 . 26 ) 

Me R.€ 

In + i 2In ^ 0^®) 

Minimising these residues by the Galerkin method, 

one gets after integration by parts, 

h. / ^ 

; dft. +Mc Re^ ^ 1^ dh 

*U 

= N? A^(Z^®^^ - 0^®) ) 


/ 0^®^'dh ^ ^ A (Z^®^^ - 0^®) ) dh 

r \ Tv .' T ^ V i il 


Me Re 


+ Re^ J wf I 0^®) dh - i2 CRe^ / N? r Ji^®^ dh 
Q i u 0 1 n 


> N? I„ 5!(®) 


1 n 


h 


(2.27) 

.hatching the above equations, it is seen that compatibility of 
Z and 0 at nodes and completeness of Z, Z* , 0 and 0^ is 
required. A linear polynomial i.e. element for Z and 0 
satisfies the above requirements, but it lea ds to i ll -conditioni ng 
[16 ] . So a cubic polynomial i.e, element .is chosen for 
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Z and 0 both, V.’ith this, fewer elements need be used for 
same accuracy. Since = N^, the superscripts are dropped 
nov; on\;ards . 

Substituting for and 0^®^ and arranging Eqns. (2.27) 

in matrix form, one gets 

h , , (ne) ^ , , (ne) 

J -n [N JciM0] 

° ° . h 

o b .. (ne) .(e). 

+ Me Re^ / A {N? [NJ dh {Z} = {NJ A (Z'^ ^ - 0*^ 0 


h (ne) ^ (^®) 

Me Re^ (n'j [N' Joh 10] - / IN} [N' J dh tZ} 

• 0 0 

h (ne) /. b •: (ne) 

+ ; A ^K] [_N Jdh ^0} + Me Re^ / ^N] [N J dh ^0} 

o . ° 

( ne ) p / I b 

- 2i Me Re^ ; I„ {N] [N j dh {i} = {N} Me Re In 0 

0 * '0 


.... (2.28) 

These equation can be arranged in. the form 

2 n(e) rn 1 *1 r/Vi C [O ] [^1 12} ^ ^ 

McRe'^[Mi^]^ ^ [0 ] fZ] ^ _j_ 

col KoReVi]'®’ [0]a.IcRe^[MI]'®j 


i ] 


71 (ne) 


-m 




(ne) 


^FBS} 


(ne) 

.(ne) 


( 2 . 29 ), 


# « « • 



where 


XI 

[ {n‘} [N' J dh 

o 

{n‘} [N J dh 


r ^ I ' I 

[KI] = /I lN}[N'jdii 

o 

= / A {Ni [Nj dh 

o 

[ mi/®) = I„ iNl LNj dh 

A„(z(®)'. 0 ) 1 

fPBS3^’^®^=J 0 

Aj,(z(®)' - 0 ) 1 



{FBI13 


(ne) f-McRe^ In 


0 

McRe^ 


.... (2.30) 


‘k 


(q\ 

It may be noted that here all matrices [ M]^ 

[G]^®^ and [K]^®^ are symmetric, 

UquEtion (2.2A) is the governing equation vihen shear 
deformation is negligible, -i^is is a single fourth order 
differential equation. Its PEM equations are obtained in the 
routine way and are 



(n®).io[o](e) [r^®) jq j(ne), j rbj (ne) 

.... (2.31) 


Where , 


i®) ^ r , 9 ^ . 

[M] = ; {nIlNJ dh + Re"^ ; JN'J [N'Jdh 
o o 


is the elemental mass matrix, 


[G 




aRe*^ /[N^HN'Jdh 


is the elemental gyroscopic matrix, 

(e) h 

[K] =» / [N"Jdh 

o 

is the elemental stiffness matrix, 

(ne) (ne) 

{1^ “ ?25 

is ttie elemental nodal displacement vector 


ana 


-in5 


+(N‘j Z 


(e)" 


■[N| Re 


2 V' 


th 


-.[N^2Re^ 


is the elemental noaal force vector. 

It may be noted -ttiat here C..| element is needed, and 
all matrices are symmetric, 

2,2 DUAL SHAFT s 

Non-^circular shafts or circular shafts With key ways 
and slots will have different stiffnesses in different 
directions, dtiffness in any direction can be deter mined, once 
stiffnesses in two principal directions are luiown, T^us if 
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ana are -ttie principal moments of inertia in principal 
directions I and II, then moments of inertia in X and Y directions 
are, Fig. (2. 3 a), [17 ] . 


XX 


^ ^ )C03 2*? 






^yy “ ^-T-> - 2* 


.... (2.32) 


xy 


(--■ " ■g-— ) sin 2!? 


Such shafts are referred as dual shafts in literature. 

For such shafts the moment curvature relations are [18 ] 

= - [EI^ 0^ + El^y ®y ^ 

.... (2.33) 

= [Eiyy e; + EI^ ] 


Shaft is rotating with angular velocity Cq. X and y 
are the fixed axes i.e. they do not change, -f*t any time the 
principal axes make an angle o^t withX,y-axes Fig, 2, 3(b). 

Using the Eqns . (2.32) (^=- C^t) and (2.33) and 
following the procedure of Appendix II, the equations of 
motion for uniform dual shaft become^ 

For bending in XS plane : 


00 s 2 B(i!^ sin 2 O^t) 


vj GA (x^ - 0y)’ - mx s 0 

I, -I 




.... (2.3A) 


+ k'oa(x'. Sy). Ijj 6y - Ip Qp e* " 0 
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1 

(a) 



I ^ (b) 


Fig:2.3 (a) Position of frame(XYS) relative to (I IIS) 
(b) Position of frame (IHS) relative to (XYS) 



2^ 


For bending in YS plane i 


K^GA (y* - © )' - my = 0 


E{h^ + 51^2 


I —I 

cos 2 Q^t)©J^-E(-l^sin 2 O^t) ©{J 

’Z 


.... (2.35) 


+ K GA (y'-e^) . Ijj + Ip 0^ 6y = 0 

Multiplying Eqn, (2,35) with i =« ^"-1 and combining with 

Eqn, (2.54), using the notations z = x+iy f Q = 0^+ i©^ and 
© = ©y - i©^, the equations of motion become 

k'ga ( z' - ©)' - mz = 0 

E(fl^) e". E (il^) cos 2 8" 

-iE ( — Ig-^) sin 2 Q^t + k‘gA (z* -©) 


.... (2.36) 




Combining second and third terms of second equation, one gets 


K^GA (z* - ©)‘ - mz » O 


E(^l^)©“- E(^^J^) 5'e^^ + K^Gii (z’- ©) 


.... (2.37) 


- Ijj © + ilp Q^© « 0 

-i 0 t 

Multiplying Eqn, (2.37) by e 


KlGA 


/ i -v' ~i D t ** -i 0 t n 
(z - 6) e 0 - mz e o = 0 


I-+I 


I. -I 


'’"'"^)e"e^ °o'^ -E(-l^)5 '' +k'ga(z'- ©)e"^°o‘^ 


, . . . * (2.3b) 


-Ij3 e e-i “o^ + ilp “p ® - 0 



25 


Equations (2,58) have ooefflcients which are functions of 
time. This time dependence can be elimin^’ted by defining 


5 = 26 


-i 


0 t 
0 


^ s Oe 


-i 




^ ss S 




(2.39) 


Thus Eqns, (2.38) become 

k‘gA (s' - mV - i2m 0^1 + mC^S = 0 

E(is^) k'ga 

.... (2. to) 

-Id - I2Id Oo^ ^ Id “o 

+ ilj, - Ip 0^ ?! = 0 

TaKinc Ip, as done in [A ], the above equations become 

k'gA(S* - mS -i2ra Q^S+ra O^S = 0 

5!'L E(il^) ^VoACs'-ili) .... (2.41) 

- Id «0 ° 

These are differential equations in complex variables 

. S=‘So''' i^o» ^ ^ (2.42) 

w.iere , h are the deflections along I and II respectively 
and the bejiding slopes in IS plane and IIS 

plane -respectively . 
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Seprir.itin^ real and imaginary parts, one gets 

k'g-C-sJ,- o*^°oV"'°o\ = ° 

Eig e" +K'GA(^'^-e„)-ijje„.i^ e„ = o 

K'o>i(’i^- )' -jnh^-2m * 0 

EI^ e^' *k'ga(»i;. Sj) -1^9^ -Ij, d| 0^ ■= 0 


.... (2.44) 


Equations (2,43) correspond to IS plane and 
Eqns. (2.44) correspond to IIS plane, 

I 

It may be noted that Eqns, (2,43) and (2,44) differ from those 
given by Dimentberg [4] . It is because of the definitions 
of I^ and Igfvhich differ from Dimentberg's definition. 

By observing Eqns. (2,39) i it may also be noted that, 
the natural frequencies obtained using these conversions, 
will be the natural frequencies (>/ ) in I II S frame , 

To obtain the natural frequencies (\) in XYS frame, the 
follov/lng relation is to be used. 

\ = \^ + o • . , . (2,45) 

o 

Two second order equations of Eqn, (2,43) can be 
combined into a fourth order equation and same can be done 
for Eqn, (2,44), These two fourth order equations are same as 
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•Lhooc given by Dimentberg [ ^ ] t v/ere oblrilned in a direct 

way \iiitl.ouL deriving the second order equations. These second 
order equations are of great advantage in FEn formulation, 
as all boundary conditions (geometrical and natural) can be 
incorporated exactly, V/hen shear deformations are negligible, 
these tVi/o fourth order equations become for bending in 
IS plane, 


^^2 '^ 0 + “o i to-2 Vo- °o ° 

(2.46) 

and for bending in IIS plane 

« 

El. + 0?’^!!) +in (h‘ +2 £} t - qJ n ) = 0 

To D o 00' O O O 0 O 

.... (2.47) 


2.2.1 Non-dlmensionalisation : 

The non-dimensional parameters used are 




Q 





C - 


.... (2.48) 


where k is the radius of gyration for mass moment of Inertia 
and is the radius of gyration for area moment of inertia I^ . 
V/ith these the non-dimensional form of Sons. (2.43) and 
(2.44) are 



2S 


( tl- 0^)- Me Re^ (t. -2 on- Q^%) = 0 
MoRs^ ij 1 A o 

— ^ ®T) + U ^ 677) “ ° ®r)) = 0 


..., (2.^9) 


and 


(r)'- /- Me Re?(i7’ + 2 0 -'^a^n) = 0 


Me 


Re^ + (n' -, 0 ^ )-C^ Me Re^ (©^ + 0^0^ ) = 0 


.... (2.50) 


? I'ii 

where Re^ = ■■ ■ « ■ , is the rotatory inertia coefficient and, 


•1 


£ 


AT 


Me ss is the material constant. 

k'g 


Similarly Eqns, (2,46) and (2,47) become 


C C^Re^ (?,"+ 0^5,^0+C (i -2 O’i- 0^^ ) = 0 .... (2,51) 

and 

Re^ (T1^+ 0^ n^^)+(Ti + 2 Q4- Q^’^) “0 ..,, (2.52) 

2,2,2 FEM Equotions ; 

Proceeding as in section (2,1 ,2^ finite element 

equations for dual shaft are 

(e) (ne) (e) . (ne) (e) (ne) (ne) 

[K ] l<fj + [0] tq] , +[k] {q} =iFB} 


• # * • 


(a.53) 
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wliorG 



Mo Re^ M/|] 

[0] 

[0] 

[0] 


[0 ] 

^Mc Re^ [ M^] 

[0] 

Me Re^[M-,] 

[0] 


[0 ] 

[0] 

[0 ] 


_[0 ] 

[0] 

[0 ] 

jMc Re 

C 1 1 

s 
« — » 

is tlie elemental mass 

matrix , 




r- 

[0 ] 

-2 0McRe|[M^] 




[0 ] 

[0] 


[0(^)3= 

[0 ] 

2 OMoRe^pi^] 

[0 ] 

[0 ] 

[0 ] 

[0] 

[0] 

[0] 



-[0 ] 

[0] 

[0 ] 

[0]_ 



is the elemental gyroscopic matrix, 


([IC| ]-McRe^ 0*^14-,]) 




Sym, 


[Kg] 

MoRe? 

(-V 

+ (1+C^jiyicRe^ 0^) 


[0 ] 
[0 ] 


2 2 
-McRef o'^ 


[0 ] 
[0 ] 

[Kg] 





(KoRe^ [Kj 
+(1+C^McRe^ 

. fi 2 ) 


is the elemental stifinees matrix, 
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(ne) _ 


iu 

{n] Cne) 


i®£.] 


(ne) 


is the nodal displacement vector of the element, 

{N} 


{EBj 


(ne)_ 


McRe 


- IN] 4 


(e)' 


^N] 0^®^ ) 


o 

h 

o 

h 


McRe^ JN/ 0^®^ 


o 

I ,h 


is the nodal force vector of the element, 
and 

h 

fM.] = / {N| [Njdh 

U I J Q 

h 

[Ki] = J iJ'’’! [N'Jdh 
0 


.. (2.54) 


h 


[Kp] =.- / ^N^iLNjdh 

C. • Q 

f \ (® ^ r T (e ) 

Here and [k] are symmetric matrices, and [G J' 

skew -symmetric matrix. 

Similarly FEM equations for a dual shaft without 
shear ^ Sqns, (2.51) and (2 .52)^, become 


is 



(ne) 

= {FB] .... ( 2 . 55 ) 


where, 

C([Mi] +C^Ref[lC,]) 

ri.ij(e) 

[ 0 ] 

is the elemental mass matrix, 


[ 0 ] 

([M,]+CjjRe2 [i^] ) 


[<?] 


(e) 


[ 0 ] 

l_2[Ml] 

is the elemental gyroscopic matrix , 


-2C [K|] 

[ 0 ] 


[K] 


(e) 


[Kj] *CC^(C_jRef [K^] 4 ;Kl,] ) 

[ 0 ] 


[Kj]* O^CC^Ref [K^] 

- [ Ml] ) 


is the elemental stiffness matrix. 


{FB} 


( ne) 


f 

0 


+ (N<|T7(e) 


h , h 

+ {Nisi ® +CC^Re^^ ^ ® ^ ^ ^ 


0 

hih 




‘0 'o 

is the elemental nodal force vector, 




ih 
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is iliG Glemenial nodal displacement vector 
and 

h 

[ll-ll = s LllJ d^' 

o 

h 

LK^J = ; 5k'J [N'jdh .... (2,56) 

0 

[K 2 ] = / {N*}[N"jdh 

o 

For both the cases element has been used as 
discussed in section (2.1 ). Here [M and [i: are 

symmetric matrices, and [g is unsymmetric matrix. 

2.5 ISOTROPIC SHAFT V/ITH DAMPING ; 

The friction offers resistance to motion and this 
usually leads to damping of natural vibrations. The friction 
arising between rotating and stationary parts ,such as bearings 
is cilled as external friction and damping arising due to 
this is called external damping. The friction arising within 
the rotating shaft due to resistance of the particles of 
material is called internal friction and the damping due to 
this is called as internal damping . 

External damping can be incorporated by modelling the 
bearing as dash-pot model v/ith a viscous damping . (Coefficient C 




•7 *7 


and e'jiK'.tinL the shear force to Lhe dampinj, force, while 
"the boundary conditions. 

To take internal damping into account, tv/o forms of 
damping moo els are taken. In the first, it is linear 
velocity dependent viscous model, for w'hich the constitutive 
relation relating the axial stress a to the axial strain 

o 

Is [ 9 ] . 

....( 2 . 57 ) 

Whereby is viscous damping coefficient. 

The other model is frequency independent hysteretic 
damping v/hich is expressed by the loss, factor which is 
related to loss angle 7^^ [19] by 


tan 7^, = T) 

ii n 

The constitutive relation for this case is [9 ] 


.... (2.58) 


H 


n 


H 




] 

' 3 


(2.59) 


H 


Combining the two models, Eqn. (2.57) and Hqn. (2,59) 
one gets 


1 


n 


0^ !■( )e^+ (T)^+ 


H 


\ JTPn, 


)eg] (2.60) 


H 




fieferrin^ to (2,4 ), z (s,t) is the distance of 
neutral a::is of cross section, K is the speed of precession, 
C o is the spin speed. The normal strain in the differential 
area dA , at a radius r j is 


= . i = .roos(0 ^ 

C o s 


(2.61) 


v/here 6 is the perpendicular distance of dA from neutral axis 
From Fig, (2,4), moments about X- and y axes are 


2n r 

= S / [ y + r sin D t ] dA 


2ti 


4 • , • (2 , 62 ) 




Using Eqns, (2.60) and (2.61), moments from 
uqn, (2,62) become -[9 ] , 


M 


X 


1+n 


H 


=EI 


Jun^ 


H 


(QPv^ 


( Co\+ 


n 


H 




■) 


r, 1+r) 

Ji_ ) ( — ±. ) 


H 




x^' 


+ El 


-n. 


'V 



J/l 


x*‘ 

L J 


(2.63) 


/\a-a A»’'i 

Using the above Eqn. (2,63) in Eqns. (*»§■) and (#=4) of 
Appendix JI, 'the governing differential equations are 
obtained as, 



[ 


. o ) yiv ] 

/Uiig luT? 


- I x'^ 


H 


“I Q y‘‘+ n6i - 0 
p o*' 


(2.64) 


for bending in XS plane and 


1 l-T), 

El [ — ri- y'^''- ( 


-i'' -Iny’' 


Ju4 


H 


+ Ip x“+ my = 0 


for bending in YS plane. 


2.3,1 Non-dimensionalisation : 

The non-dimensional form of Eqns . (2.64), are obtained 
follov;ing the procedure of section (2.1.1). The non- 
dimensional parameteiE used are, 









3 

El~ ’XX El T* 


• • • • 


(2.65) 


Ihen.Eqns, (2.64) in the non-dimensional form, are 


(_!!!IH) xIV ^ (Jk. . 0 ) ^ ’Vn 

- 



- Re‘ 


(x" 


2 oy") 


+ X = 0 
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1+71 n 

( — :^) - (— iL + Q ) + n 


:riy 


/1^ 


V.; 


K 


2 / 


.... (2.66) 


- Re^ ( Y"- 2 QX''^) + Y = 0 


where Re‘ 


^1 


2,3.2 FEi. Equations : 

Follov/ing a procedure similar to the one explained 
in section (2,1.2) in detail, FEM equations are obtained and 
expressed in matrix form as follov;s, 

{q] + [C] +[k] ^q^ =1FB} 

(2.67) 

v/here 


[M ] 


(e)_ 


[M^] + Re^ [im 
[0 ] 


[0 ] 

[M^] +Re^ [Kj 


is the elemental mass matrix, 

7iy [K2] 2 QRe^ [kJ} 


[C ] 


(e)= 


_~2 [K^ ] 


\ [J^2] 


is the elemental damping matrix. 


[K] 


(e) 


1+n 


( H) [k^] 


"h 


) [KJ 


H 


n 


1+n 


+TJyj^Q )[ K^l C" 


'H 


) [^2] 


fTn 


H 





is the elemental stiffness matrix, 




Cne) 


.fxf 


(ne) 

(ne) 


is the nodal displacement vector of the element, 


{fb5 


(ne) 




(ne) 

(ne) 


is the nodal force vector of the element, 
where 

V, 

JTThg 




,h 


( ” +"VH “ H - 

{N} Y^®^''| 



r- 

"VN^ ' ’ . 

rh , . 

+ {N ] A 

0 


h .h 

o 




,h 


+ Re^ [ X 




[h 


o 

I .h 


] . 




(ll2i ) [ - + {n'] 


h 


^/U^ 


'H 


- (r^ * "VN ) t - 




h 




+ Re2[WY(®>' 


h 

o 

h 
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( 2 . 68 ) 


Here also element has "been used Matrix 
is symmetric and matrices and [K]^ ^ are 

skew-syunetric , 



CHnFTER- III 


results akt> discussion 

\ 

Computer pro^r’ams in FORTR.-.I'' IV have been developed 
to determine the naturol frequencies and critici-l speeds, 
for all the cases of shafts, formulated in Chapter 2. 


3.1 METHOD OF SOLUTION : 

Element matrices, [g] ^ [c ^ 

were calculated using Gauss quadrature v/ith four Gauss 
points. These matrices are assembled in the usual way and 
boundary conditions applied. 


9 


3.1.1 Isotropic shaft without damping ; 

Assuming a hormonic solution of the form 

.... (3.1) 

anci applying to the assembled form of Eon. (2.29) the 


n . A ixt 
1°-! = ® 


follo-.."in^, equotion is obtained 


(n) 


(-X^ [ M ]+A\ [g] + [K])fq^^ =0 .,..(3.2) 

For critical speeds in forward precession, the critarion is 

\ = Q .... (3.3) 

and for baclward precession it is M - C . 

Applying these critaria., ^^qn. (3.2) becomes 


[ Kl W = 

o 


( [m1+ ^ 


♦ • ♦ * 


( 5 .*) 



where minus sign is to be used for forward precession and 
plus sign for backward precession. The eigenvalues Q 
are the critical speeds . The above eigenvalue problem is 
solved using li4G library subroutine 

3,1,2 Dual shaft v/ithout damping; : 

The method of solution for determining the critical 
speeds of dual shaft is given in this section. Governing 
differential eo_uations for this case are in reference frame 
I II3 , The natural frequencies (\' ) in this frame of reference 
are related to natural frequencies (K) in frame of reference 
XYS , by Sqn, (2,45), It is to be noted here that in the 
Eqn, (2,45) is the spin speed in dimensional form, which is 
later non dimens ionalised as c , 


Noting that, the critarion for critical speeas to occure 
in forward precession is K = c , in reference frame XYS , 
the critarion in I IIS frame becomes 

= 0 (3,5) 


Considering a dual shaft without shear, i-qn, (2,55), the [k]® 

matrix can be expressed in the follov.''ing form. 

, le) (e) 2 (e) 

[ K J = [KD] - Q [md] 


• « • • 


(3.6) 



A 2 


\Jherc 


[ro] 


(e) 


and 


[ Ml; ] 


[Kgl [0] 

-to] [K2L 

0([M^] -C^Re2[li,] ) [0] 

[0] Ref [k^] 


Then -i^qn, (2,55) can be re-wribben as 




[g](®) {q|'"®^+( [KD]^®)-c 2 [MD])4,5('ieL^FBf(nO 


.... (3.7) 


-•.ssuming a solution of the form 




.... (3.8) 


and applying to assembled form of Eqn,(3,7) and noting the 
condition Eqn, (3.5), the eigen value problem becomes 


[KD] = 0^ [m .... (3.9) 


in v/hich, the eigen values represent critical speeds in 
forv/ard precession and Eqn. (3.9) is solved using M 
library subroutine F02AEF, 

For critical speeds in back\‘rard precession, the critarion 
in Xys frame of reference is 


X *= — 0 


* # • « 


(3.10) 



A3 


Then, from hqn. (2.A5), the criterion in I IIo frame 
becomes 

- 2 0 .... ( 3 . 11 ) 

Assuming solution of the form given by 3qn. (3.8) and 
spi'lying to the assembled form of Son, (3,7), the eigen value 
problem, to be solved becomes 

[ KD]{q^ = 0 2 + i2 [g ] + [mD ] ) .... (3.12) 

where the eigen values represent the critical speeds in 
backv.'ard precession. Nov; as one matrix is complex, the above 
eigen value: problem is solved using the MAG library 
subroutine F02GJF, 

The critical speeds of a duel shaft with shear 
deformation effect are also obtained using a similar 
procedure , 

3.1,3 Isotropic shaft with damping : 

The nature 1 frequencies for this case are obtained 
solving Eon. (2.67), for a range of spin speeds, using the 
method of solution giveb below. 

The Eqn. (2,67), after assembling and applying 
boundary conditions becomes 

[ m] c] +[ k] = 


0 



whicii c:.n be re-written as [20 ] 


[0] [M] 

p] [4 

where = 

i'*ssumin:, a solution of the form 

(3.1A) 

v/here "K — Xj^ + fs complex, Eqn» (3»12) becomes an 
ei^isn value problem of the form 

K •••• 

The above ei^en value problem is solved using NAG library 
subroutine F02 aFF. The eigen value? appear as complex conjugate 
pairs, with positive or negative real parts. The positive 
real part indicates an instaoility , 

3.2 COlNER'^SliCZ STUDY : 

To study the convergence, results (critical 

speeds) v/ere obtained for a steel shaft of uniform circular 
section, using 4,6 and 8 elements, by solving Eqn. (3.4 ), for 


[ 0 ] 


[I] 


- -[K]"”’ [C] 


[M] [0l 

[0] [K] 


= 0 .... (3.13) 
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hinged -hin^s^j fixed -fixed and caitilever bounJary conditions. 

The volue of material constant (he = — 7 — ) for such a steel 

k' G 

shaft is 2 , 9 . ^he value of rotatory inertia co-efficient 
(Re = used is 0,05, which represents a shaft having a 

length, five times its diameter. 


The I’esults are presented in Table. 1 .For hingeci -hinged, 
boundary conditions the results are compared v.’ith the closed 
form solution given by Dimentberg [ 4 ] , For other boundary 
conditions, comparison is done \/ith the results obtained by 
God [16 ] using transfer^ FEM with a single lourth order 
governing differential equation, 

.iith an increase in number of elements used, the 
results converge from above and the error is larger in higher 
modes, as is to be expected by FSM, The critical speeds 
obtained using 8 elements were found to give accurate results 
for all the boundary condition for the number of moaes studied 
Hence 8 elements were used for all the cases studied in 
this v/orl:. 



Tabic 1 : 


Conver^^cnce ^tudy of Critical Speeds 
(kc = 0.05, He = 2.9 ) 

Hinged -Hirioed Shift 


Mode 


T 


Ho» of elements 


Exact [ 4 ] 


1 


9.6587 

9.6383 

— a- 

9.6383 

9.6383 

2 


55.8967 

35.8793 

35.8774 

35.8768 

3 


72.3132 

72.1620 

72.1429 

72.1370 

4 

113 .6066 

112.3852 

112.2958 

112.2660 

5 

156,5399 

153.5009 

153.2122 

153.1123 




Fixed -Fixed 

Shaft 


Mode 

1 

1 No. of elements 


I 

r:. 


8 

La r i6i 

1 


19.3408 

19.3379 

19.3377 

21 .660 

2 


47.0807 

47 .0430 

47.0367 

54.682 

3 


81 .5682 

81 .3770 

81 .3501 

94.827 

4 

120.4545 

118.9659 

116.8833 

136.970 

5 

161 .4977 

158.0973 

157.8044 

- 




Cantilever 

Shaft 


Mode 

"1 

jsio, of elements 

i 

- . ... -J 

L g 1 

1 5“ 

6 

8 I 

j" o L j 

1 


3.4770 

3.4770 

3,4770 

3.477 

2 


20.4524 

20.4492 

20.4490 

20.441 

3 


52.0531 

51 .9979 

51 .9915 

51 .902 

4 


90.4899 

90.2031 

90.1609 

89 .801 


1 30 .9754 


130.8163 


5 132.9002 
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3 . 3 SHAFT ..IT'-iOUT £> IIFIIIG ; 

-'.5,1 Jl-Cect of gyroscopic action and rotitpry Inertia 


'J-yroscopic action term is the [G ] matrix in 
^un, (2,29) . Neglecting xhis term, the equations become os 
those of Timoshenko beam. Results ( critical speeds ) for 
forv/ard precession for various shafts are compared witn those 
of Timoshenko beams for Re values of 0,0125 and 0,05 in 
Table 2. 


iiesalts for hinged -hinged and cantilever shaft match 
well v.'ith the exact results of [4 ] and [ 15 ] » but there is 
large difference for fixed-fixed shaft. This can be explained 
by notii\ th. t Goel’s [16] results are based on single governing 
fourth order differential equation, whereas present results are 
based on tv/o second order differential equations. In Goel’s 06 ] 
approach one con satisfy only one geometrical boundary condi- 
tion (3BC) Z = 0 exactly. Other geometrical boundary condition 
of slope and two natural boundary (r©C) conditions on shear 
and bending moment c^in be satisfied only in an approximate way 
whereas in the present approach all boundary conditions 
(geometrical and natural) can be satisfied exactly?, in case of 
hinged -hinged shaft,Goelfe approach satisfies two GBC exactly 
and two MBC approximately . .And for fixed -fixed shaft it 
satisfies two GBC exactly and other two GBC approximately. 

It is seen that approximate apr)lication of two GBC in case of 
fixed -fixed shafts leads to a large error. 



lable 2 : btudy of gyroscopic and rot-tory inertia 
effects on critical speed's i 

(He =2.9) 

Kinged -hinred 


liode ■ 

L He = 0.0125 

Re= 0, 

.05 




iiimoshenko 

f ^ 

Q Shaft 


Tiposhenkoi 

5 Shaft 


j Desrn 

ilGoel 1 

L161 

Present 

bear. 1 

iGoel 1 

■a 

Present 

1 

9.84 

9.85 

9.86 

9.43 

9.64 

9.64 

2 

39.02 

39.26 

39.25 

33.80 

35.89 

35.88 

3 

86.58 

87.77 

87.70 

6C .47 

72.22 

72.14 

4 

151 .14 

154.81 

154.43 

103.11 

112.66 

112.30 

5 

231 .11 

M* 


238.43 

I4l .54 

- 


153.21 



Fixed 

-fixed 





T*] on *1 

1 Re 

= 0.0125 


i Re = 

0.05 




jTimoshenko 

j 

Shaft 


i 

1 beam 

■*Bggwi«aii 

Present 


ea 

Present 

1 

22.11 

23.34 

22.15 

19.21 

21 .66 

19.34 

2 

60.02 

61 .29 

60.42 

45.04 

54.68 

47.04 

5 

115.42 

119.39 

116.97 

76.53 

94.83 

81 .35 

4. 

186.47 

195.84 

190.53 

110.85 

136.97 

118.88 

5 

271 .53 



279.92 

146.96 

M* 


157.80 


Cantilever 


\ Re = 0.0125 


Re = 

0.05 


H0C16 


Shaft 


\ Shaft 






■IRBSBHi 

mmmfm 


1 

3.51 

3.51 

3.51 

3.94 

3.48 

3.48 

2 

21 .82 

21 .93 

21 .93 

19.23 

20.44 

20.45 

3 

60.32 

61 .03 

61 .01 

47.14 

51,90 

51.99 

4 

116.10 118.62 

118.46 

80.15 

89.80 

90.16 

5 

187.83 


193.57 

115.91 

- 

130.82 














U’ 


It is seen that, the effect of gyroscoi-dc action 
is Lo inci-’eace the natural freq^uency of the system and the 
effect is significant in higher modes. This is hccauce, the 
moments arising due to the gyroscopic action, cry to 
straighten tne shaft, and thus account for an increase in 
effective stiffness, which increases tine natural frequencies 
of the system. The gyroscopic moments increase with increase in 
curvature, .*s the curvature is larger in higher modes, the 
effect is of appreciable magnitude in higher modes . The 
effect is seen to increase with an increase in 'Re. A larger Re 
implies a thicher shaft (i.e, a smaller l/d ratio, where d 
is the diameter of shaft), which also leads to larger 
gyroscopic moments, and hence increased effect. 


It is also seen that critical speeds decrease with 
incro:se in Re value, and this effect is significant in 
higher modes, as expected, 

3.5.2 Rffect of shear deformation ; 

Critical speeds in forward precession are obtained 
by solving Eqns . (2,31), and are given Table 3 for various, 
shafts, for two values of Re. Results with shear are also given 

in this table , 

It is clear from Table 3, that neglecting shear 
deforma.tion effect 'leads to an over-estimation of critical 
speeds and the effect is of appreciable magnitude for thicker 
shafts (i.e. higher Re), which is as expected. For Re = 0.05 



Tc.tle 3 : Study of she:^r deformation effect 
on critical speeds . 

(Me =2.9) 

Hinged -hinged 


r.odc 


nithout shear 


9.88 

39.61 

89.57 

160.53 

253.98 

518.80 

739.01 

933.83 


2:3 


vith snear 


9.86 

39.25 

87.70 

154.43 

238.44 
338.57 
453,71 
589.84 


Re = O.Op 


lout shear ijVvith shear_ 


9.99 
41 .60 
100.84 
203.81 
402.69 
1092.45 


9.64 

35.88 

72.14 

112,30 

153,21 

193.89 

234.20 

275.21 



Mode 


= 0.0 


it h out shear 


22.40 

61 .94 
122.14 
203.86 

308.94 
439.68 
592.05 
847,90 



22.15 

60.42 

116.97 

190.53 

279.92 

383.80 

500.15 

539.98 


= 0.07 


Without shear 


22.73 
65.60 
139.66 
266.20 
521 .53 
1926.30 


19.56 
47.04 
81 .35 
118.86 
157.80 
19'7.16 
236.59 


277.08 







'.T^'iRAL library 


^^cc. No 


KANfUR 

(jT2^ 
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Taole ^ ••••( contd «) Cantilever 


] 

Mode 


0T25 

1 Re 


■■oTo5 



:[ v'ith shear 

1 Without shear 

T 

With shear 



$ 

± 


1 

3.52 

3.51 

3,54 


3.^8 

2 

22,09 

21 .93 

22.99 


20.45 

3 

62.11 

61 .01 

68.72 


51 .99 

4 

122.55 

118,46 

150.87 


90.16 

5 

204.72 

193.57 

303.32 


130.82 

6 

310.38 

285,20 

679.33 


171 ,84 

7 

4A1 ,62 

392.04 

- 


212,56 

8 

593,75 

511,67 

M 


252.90 
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only six critical speeds are obtained in forv/ard precession 
for all the tliree shafts. This confirms the v/ork of 
Dimentberg [4] , who showed by close-form solution that the 

number of critical speeds (n) in forward precession is 
limited by (for hinged -hinged supports). 

n < (3.16) 

Observing the results for fixed -fixed and cantilever 
support conditions, it seems that above relation holds good 
for these shafts also. 

3.3,3 Multi -span shafts; 

Three shafts with two, three and four spans used 
by Goel [ 1 6] are studied and are shown in Table 4 , All 
supports are hinged. Number of elements used in each span are 
shovm in brackets. Critical speeds are presented in Table 4. 
Also presented are the results of rotor given by Gael 
using transfer FEM. 

With increase in number of supports, the stiffness of 
the system increases significantly . Hence, a significant 
increase in critical speeds is to be expected, and it is clear 

from Table 4, 

Except the results for first mode for three and four 
span shafts, the present results are comparable with those of 
Goel [16]** two results for the first mode are clearly 

in error as they violate the fact that increase in number of 



Tabic A : Cri Lic?l speeds for multi-srjan rotors 

( Re = 0.05 , He = 2.‘?) 
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Note : 1 ) The numbers on top are the ratio of span length/ 
total length 

2) The numbers between two supnorts^are the number of 

elements considered for that r.span, 

• * 


Two -span rotor; 



Mod e 

Kotor 


Present 

Goel [16] 

1 

23.98 

24.95 

2 

69.60 

70.62 

3 

93.87 

104.22 

A 

128.92 

138.36 


Three- 

-scan rotor : 

.. 0.4 ^ * 

^ (A) C2) I [zy « 



• Mode 

Kotor 


Present 

Goel [ 16 ] 

1 

60.10 

23.88 

2 

89.88 

86.40 

5 

113.42 

111.91 

4 

160.22 

142.26 
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Four-sijan rotor : 


o.;tg o-a 0 3 Q 

flc2) 12> A (2)1^ 


I'lOde 


Rotor 

Present 

Goel“ 

1 

100, 2S 

31.69 

2 

121.1 2 

112.63 

3 

128.78 

133.25 

4 

180.64 

162.73 
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supports increases the stiffness significantly, I-Iorever these 
types of inaccuracies are quite common in transfer FEM. 
Otherwise results differ, as Goel [16] used a fourth order 
differential equation and present work uses two second order 
differential equations . 


3,4 EFFECT OF SHAFT MASS AM3 DISK FLEXIBILITY ; 

As said earlier, most of the v/ork in rotor dynamics 
neglects the msiss of the shaft and the flexibility of the disk 
Effect of these assumptions is studied in this section. 


To study the effect of shaft mass, a steel shaft of 
1000 mm length and 40 ram diameter with hinged -hinged supports 
is chosen. This shaft has a rigid disk of 30 mm thickness and 
90 ram diameter placed at a distance of 250 mm from the left 
end. For shaft, rotatory and gyroscopic effects are taken into 
consideration, i.e. governing equations are (2.31). Shaft is 
divided into eight equal length elements. Disk is taken to be 
rigid, but has rotatory and gyroscopic effects. Thus for the 
disk at node 3, one he s [^1 



• • # * 


(3.17) 
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wl'iere is the mass of disk, 1. and I are the diametral and 
a » a p 

polar mass moment of inertia of disk, and is the spin 
speed of shaft. Equation (2,31) and Eqn, (3,17) were 

assembled together and solved for critical speeds for a range 

% 

of shaft mass. Range of shaft mass is obtained by successively 
decreasing the density of shaft and maintaining all other 
pai’arneters (length, diameter and E) constant so as not to 
alter the stiffness of shaft. Critical speeds for first five 
modes as a function of ratio of shaft mass to disk mass are 
given in Table 5. 

The critical speeds for massless shaft are obtained 
by reducing the order of assembled stiffness matrix to that of 
mass matrix of disk i.e. 2x2, using static condensation [27], 
For static condensation, see Appendix III, For for\:ard 
precession .only one critical speed exists and is shovm in 
Table 5. This is same as the close form solution given by 

Tse [22], 

It is seen for this shaft and disk, in the first mode 
there is an error of 15S */. for a shaft to disk mass ratio of 
6,5e4 and 30 '/.for a shaft to disk mass ratio of 0.844, 
Moreover massless shaft assumption gives only one mode in 
forv/ard precession. Thus one should not neglect the shaft 
mass, as error can be considerable even in first mode when 
shaft and disk mass are comparable. 



Table 5 : Effect of shaft mass on critical speeds 
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To study the effect of flexibility of disks, two 
shafts, one with single disk and other with tv/o disks, shown in 
Table 6 are studied , Disk is treated as a standard shaft 
element. For the first shaft there are eight elements when disk 
is taken as rigid and nine elements when disk is flexible. 

For second shaft with two disks, there are eight elements when 
disks are taken as rigid and ten elements when disks are taken 
as flexible. 

The results are presented in Table 6, The critical 
speeds obtained assuming the disks to be flexible are hi^er 
than those obtained considering the disks to be rigid, as 
expected. In the first mode the difference in results is 5,2 */, 
for shaft with one disk and 8,1 */, for shaft with two disks. 

3.5 DUAL SHAFT : 

A steel dual shaft with follov/ing parameters has been 
chosen for this study, 

C« ^ = 1.5 , Re.= 0.05, C. (= -%) = 1.0 
^2 ^ k^ 

Critical speeds for this dual shaft with no shear deform^-tion 
are obt>iiined from Eqns, (2.55), and with shear deformation 
from Eqns. (2.55). For without shear case, results for forward 
and backward precession are obtained and are given in Table 7, 
With shear case results, for forward precession are obtained 
and are also given in the same table , 




Mode 

One 

disk 

Two 

disks 

Rigid 

Flexible 

Rigid 

Flexible 

1 

473.1 

488.6 

420.9 

457.8 

2 

1806.9 

1910.4 

1800.5 

1911 .1 

3 

4392.1 

4498,3 

3883.8 

4231 .7 

4 

8326 .1 

8508.6 

8304, 6 

8470.8 

5 

12357.2 

12824.8 

11506,6 

12392.5 
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For forward precession, close forir solutions for 

hinged -hini’.ed shafts are avail^^ble [4 land are given in 

Table 7, Present FSM results match v/ell with these. It is seen 

that dual stiffness has doubled the number of critical speeds 

v/hen compared to the corresponding isotropic shaft. Further , 

iis cm be seen from the eigen v"lue problem Eqn, (3.9), f^r 

forv;ord mode of precession, the vibration in the two principal 

planes (13 and IIS) become uncoupled. Since the elas'tic 

properties are different in the two principal planes, there 

will be a critical speed critical--speed corresponding to the 

elastic property in each plane, loading to doubling ox crioicel 

speeds in comparison to isotropic shaft, 

« 

The critical speeds in b^^ckward mode of precession are 
also presented in Table 7. Also presented are, the critical 
speeds from [4], for hinged-hinged support conditions. By 
observing the eigen value problem, Son. (3.12), it is seen 
that the vibrations in two planes are coupled by the [g] . 

matrix. Multiple criticca speeds are seen to exist in 
backward mode, in contrast to foim'/ard mode. 


ISOTROPIC SHAFT '..'ITH DIFFl n Slg SUPPORT C0:'!iJ>ITI0 £ fo 
IIJ TWO PLAIila'^; 


In this section, the results are obtained for an 


isotropic shaft with different 
Solution is obtained from equa 


sui' ort conditions in tv/o plane o . 
tions of dual shait, iiiqn, (2,53), 
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Table 7 : Critical speeds of a duel shaft. 
( 8-elements) 


'i.ithout 
' (C= 


shear deformation (forv;ard pression) 
^1 

^ = 1 0.05, C^= 1.0) 



Kinged -hinged 

Fixed -fixed 


Mode 

Present 

Exact ^ 

uanxix ever 

1 

8.16 

8.16 

18.57 

2. 89 

2 

9.99 

9.99 

22.75 

3.54 

3 

33.96 

33.95 

53.56 

18.77 

A 

4l .60 

41 .59 

65.60 

22.99 

5 

82.34 

82.30 

114, 03 

56,11 

6 

100.84 

100.71 

139.66 

68.72 

7 

166.41 

165.74 

217,36 

123.16 

8 

203.81 

202.99 

266.20 

150.67 


With s 

hear deformation 

(forv;ard precession) 


(Mc= 2,9» Y” = 1 .5 

, Re<j — 0.05, 

= 1.0) 

Mode 

Iliriged -hinged 

Present Sxact 

Fixed -fixed 

G:antilever 

1 

7.96 

7.88 

16.57 

2.86 

2 

9.6^ 

9.54 

19.34 

3.46 

3 

30.61 

29.44 

41 .87 

17.31 

4 

35.88 

34.55 

47.04 

20.45 

5 

64.07 

69.82 

74.61 

40.70 

6 

72.14 

67.74 

81 .35 

52.00 

,7 

103.22 

94,33 

111 .50 

82.18 

8 

112,30 

103.79 

118.88 

90.20 


Table 7. ( . . .Contd, , j V/ithout shear deform-<tion (3ac’.r-'ard precession' 
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uGiHij ~ ^ • This sh'ift h-s GUi-s'Orts in 13 plane 

and fixed supports in Ila plane. Using these bound c-rp,? 
conditions the forward critical speeds are obtained and are 
given in Table £. Also given are the critical speeds, for 
similar fixing conditions an two planes i.e. similar fixing 
conditions of hinged and fixed supports . 


It is noted that critical speeds of different 
support conditions, HF-HF, can be obtained by combining the 
critical speeds of shafts of isotropic support conditions of 
IIH-HH and FF-FF. Thus, there- are tv;o critical speeds in 
comparison to a shaft having similar fixing conditions in the 
two planes of deflections. Critical speeds due to vibrations in 
each plane correspond to the supp-ort conditions in ChB 
respective planes. H6nce it may be said that an isotropic 
shaft having dissimilar support conditions, behaves in a 
manner similar to a dual shaft having similar support 
conditions, with regards to critical speeds. 



Sk 


Table S : Forv.'ard critical speeds ^vith 

dissirniljr sup r-ort conditions . 

(Me = 2.9, Re= 0.05, C elements) 



H: H'nged support ; F : Fixed supports 


‘lode 

HH-HH 

HF-HF 

FF-FF 

1 

9.64 

9.64 

19.34 

2 

35.88 

19.34 

47,64 

3 

72.14 

35,88 

81.35 

4 

112.30 

47.04 

118.88 

5 

153.21 

72.14 

157.80 

6 

193.89 

81.35 

197.16 

7 

234.20 

112.30 

236.59 

e 

275.21 

116.86 

277,08 
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^•7 I'/OTilOlIC 5K..FT \;IT:] LAilPIw: 


io study the efrect of internol viscous and 

hysLeretic damping, the Eqns, (2.67) are solved. A steel 

clioxt v;ith iie = 0.02 (i.e, ^ = 12,5) studied by {^9} i ^■''ith 

hinged *-hiuj;,ed support conditions is chosen. The bearings 

arc asLined to have uncoupled stiffness and damping 

charoctcristics i.e. i: = K = O and C = C =0. 

bey yx yx 

'fho conij^lex frequencies are obtained in tlie form of 


X = (3,18) 

The nondimens ional frequencies are plotted in 
i’’i^r. . (5,1) to (3,5), for a speed range of 0 to 17 (i.e, 0 
bo 120J0ri'm), and for various combin'-tions of, bearing 
flexibility and bearing damping charecteristics , 

The curves labelled ** and * B represent forward 
and backward precessional modes respectively . In all the 
crises .studied, the results are compared vrith the results of 
[ 9 ] . The observations made in each case are listed 

below. 


Case 1 : Viscous damping, isotropic undamped bearings. 


Bearings are undamped and have isotropic support 
conditions of 

Kxx" Kyy “ 

Considering the internal damping to be viscous, with 
» 0,01621 (p, 0002s), natural frequencies obtained are 
given in Fig, (3.1) as function of speed o . The 
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Fig* 31 Natural frequency Aj vs. spin speed A 
for case 1. 
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Fig.3-2 Natural frequency vs. spin speed si 
for case 2 
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0 4 8 12 16 20 24 


ja 

Fig.3.5 Natural frequency Aj vs spin speed 
for case 5 
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XirEi criticnl speed occurs at 6 .kh (49er> ri'in) and second 
crit.ice.l speed occurs at 13.55 (10500 rpm) in'forv-mrd mode 
of precision, ‘fhese values are in good agreement v.’ith those 
given by [9] , which are 6.40 (4950 rpm) and 13.55, 

The first mode v/as found to become unstable at the first 
forward critical speed and the second mode at second 
forward critical speed. These modes stay unstable upto the 
stalled speed range (15.5). 

Case 2 : Viscous damping, isotropic damped bearings. 

Tlie system studied in case 1 is analysed incorpor -ting 

isobrojlc bearing dampings with = 0.269 

(i.e, 1.75 X 10^ where C , C represent the bearing 

ui j y 

darnpin:^ coefficients in non-dimensional form in X and J 
directions respectively. The results are plotted in Fig. (3. 2). 
It is seen that, the first moae in nov: stable upto speed 
of 11 . 6(9000 rpm) and the second mode is stable for the 
entire spin speed range studied (15,5), The results match 
very well with [ 9 }- 

Cane 3 : Viscous dampixig, anisotropic updamped bearings. 

Shaft has viscous damping “ 0.0162, Bearing are 
undamped and support conditions are anisotropic i.e. 
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Iho roL.uJlts are sliown in FiZ» (3,3). It is seen th?t the 
sysbciTi becomes stable for the entire s. in speed range 
studied. Forward and b'^ckv/ard preceesionel speeds are wide 
ar.art now. Results match very well i/ith [g ] . 

C ise A : Hysteretic damping, isotropic undamped bearing , 


The effect of hysteretic dampin^ is studied in this 
c.se. Be':rin332re undamped and hsve isotropic support 
conditions 



33.14 


The value of loss factor used is 0.0002. The results are 
plotted in Fig. (3.4). It is seen that both forv/ard 
processional modes become unstable and both backward mooes 
are stable for the entire spin speed range studied. As 
be Core forv/ard and backward precessional speeds are quite 
close to each other in the modes studied, Results match 


very well with [9 ] . 


Case 5 ; Hysteretic damping, isotropic damped bearings . 

Isotropic damped bearings v/ith following properties 

are used. 


K = 

K , = 33.14 

XX 

yy 

C = 

c = 0.269 

XX 

yy 
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(/■aluo of loss factor used is 0.G002. 

Jho tl salts are shown in ^ig, It is seen xhat bot'n the 

modes are stable for the entire speed ran^ studied . The 

results match well with Fq T ■ j 

ly J • Forward and backward preeessional speeds 

ore (^uite close to each other. 


Internal viscous damping is seen to have destabilising 
effect beyond the first critical speed in forward mode. The 
stability of the system is improved by providing damped 
bG;'.rin^_,s . Stability is improved further by the provision of 
bev rings having anisotropic stiffness properties , Internal 
hystoretic damping destabilises all the forward modes over the 
entire speed range. The provision of damped bearings is again 
seen to improve the stability of the system. 

3.7.1 Effect of bearing^ stiffness on stability : 

For the shaft studied in section ‘'3, 7), case 1, with the 
isotropic bearings, the instability was found to occurs at the 
first critical speed. This case was further extended to analyse 
the effect of bearing stiffness on the instability threshold. 
InsUability threshold speeds are obtained for a range of values 
of isotropic bearing stiffness. The results are given in 
Fig. (3.6). It is seen that, the inctobility threshold increases 
with increase in bearing stiffness significantly in the begihing. 

3#7.2 Effect of bearing, damping on st ability : 

The shaft studied in Section (3.7.1), is analysed to 
study the effect of variation in bearing damping properties at 
a few chosen bearing stiffness values, on the instability 




Fig, 3.6 Variation In instability threshold 
with bearing stiffness 




unstable 


ij 
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throGhold speed, for the first noc^e. The results obf'iined 
plotted in Fig, (3.7). It is seen that v/ith an increase 
in hearing damping coefficient the instability threshold 
increases linearly. The instability threshold is seen to be 
decreasing with an increase in stiffness of bearings, when 
bcarin^ damping is present. Thus a lo’.:cr value of bs ring 
stiffness, 30.0, and higher value of bearing damping, 


0.35 (2.28 X 10^ are seen to give the highest stable 

speed range in the present case. But as can be seen from the 
Fig, (3.7^ for a lower value of bearing damping coefficient, 
h bearing v;ith higher stiffness gives better stability. 



CONCLUSIOKTS 


The following conclusions emerge out of the present 

FEM Is easily applicable to more familiar differential 
enuations of rotor dynamics. It is simple and yields 
reliable results. Satisfaction of boundary conditions 
of intermediate supports (multi-span rotors) is 
s braight forward . 

V/hen shear deformation effect is considered, results 
from two simultaneous second order partial differential 
equations using cubic polynomial approximation are 
very accurate (there is no ill conditioning) , 
Geometrical and natural boundary conditions are 
satisfied exactly. Results obtained are accurate with 
fewer elements . 

neglecting shaft mass can lead to large errors. The 
mass of the shaft is taken care of in a routine v/ay 

in FEM. 

For undamped isotropic bearings, the instability 
threshold increases vath increase in bearing stiffness 
significantly in the begining. 
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Foi' lov/er v£.lue of bearin- dainpin; higher 
stiffness increases the stability', and this 
effect is reversed for higher values of bearing 
damping , 
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-■■P2i:i€)IX - 1 

DERIVATION OF KIIETIC EIUT^Y 

The shaft cross section has a roxc-tion 0 about 
X-axis, 0y about Y-axis, 0^ about S-axis, translation 
X clcHA X-axis and y along Y-axis, as mentioned in section 
(2,1). The kinetic energy due to rotation is obtained 
as follows 

Consider a frame of reference abc attached to shaft 
cross section such thet, c axis is normal to the cross section 
aL,b axas arc in the ilone of cross section and mutually 
porx'cndicular to each other and to c axis. Then the rotations 
of crocs section relaxive to XYS frame are also the rotstions 
of frame abc relative to XYS. Initially, the axes of abc 
frame are considered to be coincident ’with XYS frame. The 
finnl position of abc frrme relctive to >2rS is determined by 
the Xollov.'ing sequence of rotations, 

1) A rotation by 0^ about X-axis gives a^b^Cg 

2) A rotation by 0^ about b 2 -axis gives a^jb^c^ 

3) A rotation by about c.|-axis gives abc 

doferinu to Fig, (A-1a) the position of a 2 b 2 C 2 ' frame 
after the rotation relative to XYS frame is 



Fig. A-1 (a) Rotation by 0^^ about X axis 

(b) Rotation by ^ about axis 

(c) Rotation by 6 ^ about axis 

(d) Final position of abc frame 
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-sin©,, 


0 


sin©. 


X 


cos©. 
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(A-1.1) 


oirnil'irly the position of freme relative to 

Cob^c^, often the rotation © is, I^i". (A -1b), 

C» ^ j 




cos©., 


sin©. 


0 

1 

0 


-sin©. 


0 


cos©. 


\ n 


CA-1 .2) 


Mnr]ly, the position of frame abc relative to frame a^b^c-j , 


after the rotation ©^ is, Fig. (^--Ic) 


a 

i b 


= 


t- 


cos©.. 


-sin©„ 


0 


sin©„ 


cos©.. 


0 


0 

0 

1 




CA-1. 3) 


The "inal position of abc after the above rotations is shown 
in Fig. (A-ld) . The frame abc has the an_ular velocities 
along -a^, see Fig.(A-^, 0y along b.^ and ©^ along c. 

The components of 9^ along the directions a,b and c can be 
obtained by the following transformation using Sqns . (A-1.2) 

and (A-1 .3) 
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oimil rly the components of 0y along the directions a,b 
c :i\. obtained using (A-1 .3) as 
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The an:,ular velocity vector is alon^ 
-lO oraiisforiB' ti on is needed. Hence tlie 
along a,b and c directions, due to the 
and ©2 are obtained using Eqns. (A-1 .4) 


c directioi'i, thus 
cngjlar velocities 
rotations 9^, ©^ 
and (A~1 .5) as 


Ca = 

®x,a 

+ 


- 9y, cos©y cos9^+©y 

Slne^ 

II 

®x,b 

+ 

®y.b+®2,b = 

K ®y 

cos© 

On = 

®x r 

+ 

©„ _ + ©_ . = 

- ©„ sin© + © 


C 

X 1 c 


y j c z ^ c 

X y z 



(A-1 .6) 


..s the directions a,b and c are also the ; rincipal directions, 
tl'iC hinotlc energy due to trsnsl-tion rot-tion becomes 


o 


(A -1.7) 


Substituting Eqn, (A-1 ,6), the eiroression for kinetic 
energy becomes. 



1 

; 

0 


[ 


m(x^+y-)+ + ®y) 

+ IpC©^ sin^©y+6|- 2 Q ^©2 sin©y) ] ds 


(A-1 .8) 


Assuming Q^, 0y to be small and noti>^.g that ©^ = % (the 
spin speed of shaft), the above expression simplifies to 

Tgl t [m(ic2+y2)+ ^ j 

f Ip( - 2-0„ Sy) ]ds (A-1 .9) 
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APPEIt)D:- II 


DEUIVATION OF DIFFEREICTL'.L ECUAi'IC.^ 3Y Ill'^^TOllS I-IETHPLi 


Free body diagram of a differential element of shaft 
is shov/n in the Fig, (2.1).NeVvi:on»s law gives 

i ( V “x 

aiy = ioy; {My+dMy)^+ 

■ Ie ®y 

51'y = nay * (Qy+®y)-^y 

®x ° ^'^x ' ”x" |2. - ^ ds 

= ((-^D “o ®y '>) 

In these equations K represents bending moment 
and u represents shear forces. Suffixes have the usual 
engineering sense in the right hand coordinate system, other 
terms have been explained in Chapter 2. Ang).aar momentums 

are chovm in Fig.(A-2.) 

Neglecting higher order terms and rearranging, 

Eqns. (a-2^)and (A-^2)become 


= mds X 

( V^^x^ % + "^x 
Ip"*® “o ®x) 

= rads y 


I 
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= mx 


ix. + Q, 


(A-2.3) 


‘V 


= In + I Q Q 

D V r) “r ■ 


X 



M' - Q = 

X 


“ ®x 


Q & 
“o y 


From mechanics of solids, it is known 


\ ' k'ga (|| : » ) 
Qy = k'ga { II - 0^) 


(-- 2.4) 


(A-2.5) 


= - El 0^ f My = + El 0y (A-2.6) 

Usinc these relations, Eons. (A^) and (rx-2,4)'become 
[ k'ga (|| - 9 )]' - idx = 0 

(A-a.7) 


(El q;)' + K'GA(|f - 0y). 0, 


0=0 

X 


[K^GA (|| _ 0^)] - my = 0 (A~2.8) 

(El k'ga(|| - 0^)- InVlp Oo ®y = 0 

Equations (A -2.7) to (A-2jB)are exactly the equations 
' (2.9) and (2,10) obtained by Hamilton’s principle. 
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A 

H. 


FPEi'IDK- 


III 


STi*TIC COKEjEIEATION 


’..hen the mass of a system is assumed to be concentrated 
at one or few chosen points, the mass matrix will be diagonal. 
Further, it contains zero diagorial elements for those 
degrees of freedom, for v;hlch there is no associated mass, 
V.'hereas the element stiffness matrix formulated by FEM approach 
is fully popul:=.ted and relates all the forces to corresponding 
degrees of freedom. But for those degrees of freedom, for 
v;hich, the mass matrix contains zero element, the associated 
inertia forces are zero. Hence, it is necessary to exclude 
those degrees of freedom from stiffness matrix, for dynamic 
analysis. The process is called as static condensation 
C Cough and Penzein [ 21 ]) and the procedure is as follows. 

Consider a typical matrix equation, 





(A- 3.1) 




where gy} represents the column of degrees of freedom to be 
condensed because the corresponding inertia forces are zero . 
represents the inertia forces associated with the 

degrees of freedom 
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From Fvin, (yi_5,1) 

r*Suc] + [Kjcy] 5y} = 

(A- 3.2) 

[Kyx^ ^ [\y] ^y} =-CO} (a-3.3) 

From Ain. (A-3.3) 

ly) = - tv? ?*} 

^ubGtitutini^ this ^yj in Sqn. (A-3.2) 

( [Kxxl - [*Scy] Pyx] > = ^'■x} 

U» 3.5) 

Thus the degrees of freedom ^y} are condensed and expressed 
in terms of ' {x] and the condensed stiffness matrix [K^] 
is obtained as 

[Kol = u 'Sex] - [K^]pyyf'‘ Pvx] > 
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